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Abstract. We study the search kinetics of an immobile target by a concentration 
of randomly moving searchers. The object of the study is to optimize the probability 
of detection within the constraints of our model. The target is hidden on a one- 
dimensional lattice in the sense that searchers have no a priori information about where 
it is, and may detect it only upon encounter. The searchers perform random walks 
in discrete time n = 0, 1, 2, . . . , A^, where N is the maximal time the search process 
is allowed to run. With probability a the searchers step on a nearest-neighbour, and 
with probability (1 — a) they leave the lattice and stay off until they land back on the 
lattice at a fixed distance L away from the departure point. The random walk is thus 
intermittent. We calculate the probability P/v that the target remains undetected up 
to the maximal search time N, and seek to minimize this probability. We find that P/v 
is a non-monotonic function of a, and show that there is an optimal choice aopt{N) of 
a well within the intermittent regime, < aopt{N) < 1, whereby Pn can be orders of 
magnitude smaller compared to the "pure" random walk cases a = and a = 1. 
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1. Introduction 

Everybody searches for something: predators look for prey, prey forage, human beings 
look for better places to work or for lost keys or for partners, spies search for hidden 
secrets, and knights desperately seek the Holy Grail, whatever it may be. 

The search for a desired target may be long and uncertain; targets may be sparse, 
hidden or difficult to detect even when found, or they may have their own life-time and 
vanish before they are detected. That is why efficient search strategies appropriate to 
specific types of situations are highly desirable. The question of optimal search strategies 
has motivated a great deal of work within the last years [1-12]. 

While earlier work has focused on systematic searches in organized human activities 
such as rescue operations [1-4], more recent analyses have been devoted to random 
strategies [5-11]. In particular, it has been realized that in the case of the so-called 
non-destructive search, when the target reappears after some time at the same location, 
and when the searcher always remains within the system. Levy flights with randomly 
reoriented ballistic trajectories provide an optimal search efficiency [5-9]. On the other 
hand, following the observation of trajectories of foraging animals such as lizards or 
fish or birds characterized by two distinct types of motion - fast relocation stages non 
receptive to the target and relatively slow reactive phases when the target may be 
detected [13, 14] - another type of random search, an intermittent search, has been 
proposed. Analytical modeling of such random intermittent search strategies has been 
put forth in [11] (see also [12]). 

In this paper we discuss the search of a single target hidden on a one-dimensional 
regular lattice by a concentration of searchers performing intermittent random walks. 
In particular, at each tick of the clock each of the searchers chooses between two 
possibilities. One is to jump to a nearest-neighbouring lattice site, and the other is to 
leave the lattice and fly with a given velocity until it lands back on the lattice at a flxed 
distance L from the departure site. The term "hidden" means that the searchers can 
detect the target only upon landing directly on the target site at the end of a flight, or 
via a one-step (nearest-neighbour) walk onto that site. Intermittency has been invoked 
in other models (e.g., see the models proposed in [11, 12]) and has been observed in a 
number of ecological and other contexts (see, e.g., [13, 14]), but all the other features 
of our model are, to the best of our knowledge, new. 

Contrary to previous work, which focused on the analysis and optimization of the 
first passage time to the target from a given location, or on the number of targets 
encountered [5-11], here we study the behaviour of a different property, namely, the 
probability that a single target remains undetected up to the maximal search time N. We 
can not, of course, make this random search process certain; that is, we can not guarantee 
that the target is found (or not found) with unit probability in a finite time, but what 
we are able to show is that the non-detection probability is a non-monotonic function 
of the parameter a which determines whether a searcher takes a nearest-neighbour step 
(probability a) or a long flight [probability (1 — a)], and that it has a sharp minimum at 
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some value of a. Consequently, we show that the search efficiency can be dramatically 
enhanced, by orders of magnitude, by choosing an appropriate value of a. 

The outline of this paper is as follows. In section |21 we describe our model of 
intermittent random walks, present definitions, and delineate our main results. In 
section 01 we derive a general formula for the non-detection probability. Section 0] is 
devoted to the analysis of different properties of the intermittent random walk. In 
section El we focus on a particular exactly solvable case, namely, a Lindenberg-Shuler 
intermittent walk with steps of unit length and flights of two lattice spacings. For this 
example we discuss the specific features of the optimization procedure. Next, in sectionlHl 
we consider the general case of intermittent random walks with steps of unit length and 
flights of fixed length L, and define an optimal search strategy. Finally, in section [71 we 
conclude with a brief summary of our results and an outlook of future work. 

2. The model, basic notation and results 

Consider a one-dimensional regular lattice of unit spacing containing M + 1 sites s. This 
is the substrate. At one of the lattice sites, say at the origin s = 0, we hide an immobile 
target, and only we know that it is there. Then, we randomly place K searchers in 
the lattice under the constraint that none of them is placed at the site with the hidden 
target, that is, we distribute them randomly over the remaining M sites. The searchers 
have no knowledge of the location of the target. 




Target 



Figure 1. A sketch of a 7-step trajectory of an intermittent random walk, starting 
at s(0), witli nearest-neiglibor steps (probability a) and off-lattice flights [probability 
(1 — a)] over a distance L (here, L = 7). The arrows indicate the step direction. Note 
that this trajectory does not find the target. 

Next, we let the searchers move according to the following rule (see figure^: At 
each tick of the clock, n = 1,2,3, N {N is the maximal time the search process 
may run), each searcher selects randomly between two possibilities. With probability 
a, the searcher jumps to one of its nearest neighbouring sites, with equal probabilities 
to the right or to the left. With probability (1 — a), the searcher leaves the lattice and 
flies off-lattice with a given velocity V until it lands at a site an fixed integer distance L 
away from the departure site. The direction of the flight is chosen at random, with equal 
probabilities to right or left. The time the searcher spends off-lattice during a flight is 
T = L/V, a model parameter that can have integer values T = 1,2, ... ,L. Note that 
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this value defines the velocity V. 

As stated earlier, the term "hidden" means that a searcher can not perceive the 
target when it is off-lattice and may detect it only when it lands on the target site. 
The searcher may land there at the end of a flight or by a one-step jump. Once the 
searcher lands on the target site, one can define an elementary detection probability q. 
The target may immediately be recognized when first encountered by a searcher (g = 1, 
perfect detection), or detection may take place with a given probability < g < 1. In 
the latter case, when any of the searchers arrives at the site occupied by the target, 
there is always a possibility that the target remains undetected (imperfect detection). 
Here we focus on the perfect detection case. Some aspects of the search kinetics in 
the imperfect detection case for the one-dimensional system will be discussed briefly in 
section El An extension of the model to c?- dimensional spaces and some other subtle 
questions not pursued here will be discussed elsewhere [15]. 

Our first goal is to determine the probability that the target remains undetected 
up to time A^, the maximal time the search process is allowed to run, which merely 
depends on our patience or on experimental constraints. We show in the next section 
that in the perfect detection limit the probability P/v that no searcher has reached the 
target up to time A^ obeys the standard formula (see, e.g. [16,17] for more details) 

= exp (-p5;v) , (2.1) 

where is the expected value of the number of distinct sites visited on the substrate 
by any searcher in an A^-step walk (see, e.g., [18]). This is a crucial equation, since we 
will arrive at results for P/v via calculations of ^iv- 

We will calculate Sn explicitly for the intermittent random walk and determine its 
asymptotic behaviour analytically in the large- A^ limit. Using these results, we set out 
to show that there always exists a value a = aopt{N), < aopt{N) < 1, which maximizes 
the expected number of distinct substrate sites visited by an intermittent random walk 
up to time A^. Consequently, it minimizes the probability in equation 1)2.11) . We 
demonstrate that by choosing this optimal value of a the non-detection probability Pat 
in equation ()2.ip can be made orders of magnitude smaller than those of the "pure" cases 
a = 1 (nearest-neighbor steps only) or a = (flights only), for which asymptotically 



and 



PAr(a; = 1) = exp 



P/v (a = 0) = exp 



-p — +0 



[2.2) 



(2.3) 



The result ()2.2|) is well-known and reflects the fact that a random walk in low dimensional 
systems {d = 1 and d = 2) compactly explores the space; that is, Sn grows sublinearly 
with time A^, and each site ever visited by such a walk is most probably visited many 
times. Hence, searching a target in a low dimensional system with the help of a standard 
random walk is not very efficient, since the walker wastes a lot of time revisiting sites 
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that do not contain the target. The result (j2.2j) is the same with the replacement 
— > N/T (i.e., it is the same as that of a nearest neighbor random walk for a shorter 
time). This further decreases the probability of detection because it does not explore the 
local regions at all before making a jump to a site that has probably been visited many 
times before. It is thus not an efficient search strategy either. We find that although 
the sublinear growth of is unavoidable in any one dimensional random walk, and, 
even more specifically, the asymptotic dependence S^- oc N^^'^, for fixed L it is possible 
to optimize the search strategy via an intermittent walk with a particular choice of a 
that depends on the length of the walk (calculated explicitly later). We obtain the 
asymptotic result 

Pat (a = aopt{N)) ~ exp 

where jiL) is an increasing function of the fiight distance L. Note that the leading term 
in the exponent contains a pre-factor (LVY^'^, V = 1,2,3,...,L, which can make it 
considerably larger than the corresponding terms in equations ()2.3|) and ()2.2|) and thus 
may result in a substantial increase of the detection probability. 



-p{LV) 



1/2 



(2.4) 



3. Non-detection probability 

In this section we define the probability P/v that the target is not detected up to time A^, 
and express it using well-known properties of random motion that have been extensively 
discussed in the literature (see, e.g., [17,18] and references therein). In particular, we 
begin by introducing the notation s{n) to denote the random process that describes 
a searcher's trajectory, and P{s\s{0);n) to denote the probability that s{n), which is 
initially at site s(0), is at site s at time n. We find it useful to present the derivation 
because, contrary to most studies, we are interested in defining the number of distinct 
sites visited by s{n) only on the substrate, which is a subset of all available space. 
Indeed, most existing work concerns random walks that never leave the lattice so that 
the equality P(s|s(0); n) = 1 holds at all times. In our case this probability is not 
conserved, P(s|s(0); n) < 1, since the searchers spend some portion of their time 
off-lattice. We present this analysis in the most general case, assuming only that s{n) 
is a homogeneous process. 

Suppose there are altogether K searchers labeled k = 1,2, . . . , K, and let Sk{n) 
denote a given n-step trajectory of the fc-th searcher. Then, in the perfect detection case 
the indicator function ^{s,fc(A^)} j of the event that the target has not been detected 
up to time for (arbitrary) given trajectories of all K searchers can be written as 

K N f 1 _ n 

^(K(iv)})=nn[i-^(^^H)]' ^(^)= ^^-^^ 

fc=ln=0 I ' ^ ■ 

This indicator function has the correct behavior, namely, for this particular set of 
trajectories \1/ = if any searcher hits the target and \1/ = 1 if none do. All searchers are 
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assumed to move independently. Averaging \1/ |^{sfc(A^)} j over all trajectories and initial 
positions thus yields the probability Pn that none of the searchers has ever reached the 
site s = up to time N, 



N 



N 



s{0) 



.n=0 



I sin] 



K 



1 



1 r ^ 



1 - / sin 




-1 K 



(3.2) 

s(0) k n=0 

where -Es(o) {• • •} stands for the average over all trajectories starting at site s(0), and 
the prime on the sum indicates the exclusion of the target site s(0) = 0. Note that 
the summation extends over all non-target sites of the substrate, a one- dimensional 
lattice, and does not include off-lattice contributions. In the thermodynamic limit we 
set K, M —>■ oo keeping their ratio fixed, K/M = p = density of searchers, to obtain 

N 

(3.3) 



N 



exp 



-pY,Esio){l-l[\l-lis{n 



s{0) 



n=Q 



Since s{n) is a homogenous process, equation (jH.Hj) can be rewritten as 



P. 



N 



exp 



N 



(3.4) 



I s=l \ n=Q 

where now Eq {. . .} denotes an average over trajectories of an auxiliary process s'{n) = 
s{n) — s{0) which starts at the origin at n = 0. Note that the indicator function 



n 



N 

n=0 



l-I s'in 



shows whether the site s on the substrate has ever been 
visited by the process s'in), n = 0,1,2, N. Consequently, the sum over all substrate 
sites is just the reahzat ion- dependent number of distinct sites visited by an individual 
searcher, and its average leads directly to the crucial equation ()2.1|1 . 

For this proof to be useful in calculating Pjy we next need to calculate S^- For this 



purpose, note that the function in the exponent of ()3.3p . 



n 



N 

n=0 



li sin] 



IS 



the indicator function of the event that s(A^) has at least once visited the origin, so that 



N 



n h 



I sin] 



E'^iv(Ok(0)), 

.(0) 



(3.5) 



s(0) L n=0 

where i?Ar(0|s(0)) is the probability that an arbitrary iV-step random process s(A^) has 
visited the origin at least once. We can therefore write 



P,v = exp -p5^'i?Ar(0|s(0)) 



<0) 



Next we note that, by definition, i?7v(0|s(0)) can be formally represented as 



N 



RN{0\siO)) = Y,Fn{0\siO)), 



(3.6) 



(3.7) 



n=0 
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where F„(0|s(0)) is the probabihty that s{n) arrived at the origin for the first time on 
the n-th step, given that it started at site s(0). Applying the theory of recurrent events 
(see, e.g., [18]), the event that s{n) is at site after n steps can be decomposed into the 
n mutually exclusive events "s(n) first arrived at site after j steps, and subsequently 
returned to site in [n — j) steps." This allows us to write for the distribution function 
P(0|s(0);n) the discrete "integral equation" 

n 

P(0|s(0); n) = (5o,s(o)^n,o + ^.(0|s(0))P(0|0; n-j), (3.8) 
which implies that the generating functions 

oo 

P{s\s{Oy,z)= ^P(s|s(0);n)z", 

n=0 

oo 

F{s\s{0)-z)= 5^P„(s|s(0))z" (3.9) 

n=0 

at s = are simply related to each other, 

P(0|.(0); z) = P(0|.(0); z) - j^^y (3-10) 
It then follows that 

oo ^ 

R{0\s{Oy,z) = 5^Pjv(0|s(0))z^ = ^P(0|s(0);^), (3.11) 

n=0 

and 

Siz) = f S,z- = ± P(OKO); .) = ^ y if;;^ (3.12) 

n=0 s{0) \ \ ^ ) 

Hence, in order to determine P/v in equation ()2.H) via S{z) we only need to determine 
P(0|0;2;) and the normalization ^sP(s|0;2). We proceed with this determination in 
the next section. 



4. Properties of the intermittent random walk 

Since the searchers in our model move independently, it suffices to focus on the properties 
of the walk of an individual searcher. The probability P(s|s(0); n) of being at site s after 
n steps obeys the recurrence relation 

P{s\s{0);n) = - [P(s- l|s(0);n- 1) + P(s + l|s(0);n- 1)] 

+ [^(^ - ^I^(O); n-T) + P{s + L\s{0); n - T)] (4.1) 

for n = 0, 1, 2, . . . , A^. Jumps between nearest-neighbouring sites occur in one unit of 
time, while long-range jumps over distance L require an integer time T. Equation 1)4.11) 
therefore defines a non-Markovian process with a memory. Since the intermittent 
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random walks are homogeneous so that P(s|s(0);n) = P{s — s(0)|0;n), without loss 
of generality we henceforth set s(0) = 0. 

The Fourier-transformed generating function 

^k, z) = Y^ exp{iks)P{s\0; z) (4.2) 

s 

can easily be calculated by multiplying both sides of equation ()4.H) by exp(iA;s) and 
and summing over all substrate sites s and time n. We readily find that 

^{k, z)=\V- az cos(fc) - (1 - a)z^ cos(A;L)] . (4.3) 

Consequently, the lattice Green function (or site occupation generating function) of the 
intermittent walk is given by 

, , 1 cos(ks)dk 

TT Jo 1 ~ a;zcos(A;j — (1 — a)z^ cos[kL) 
One of the two quantities required for S{z) in equation ()3.12|) follows immediately 
from equations ()4.2|) and ()4.3p by setting k = 0, 

P{0\s; z) = ^{k = 0,z) = {l-az-{l- a)z^)~^ . (4.5) 

s 

This leaves only the evaluation of P(0|0; z) for the calculation of S{z) and consequently 
of Sn and P/v [cf. equation ()2.1|) ]. We shall return to this calculation and to our goal of 
maximizing the detection probability of the target presently, but first we briefiy digress 
to obtain some associated interesting properties of the underlying distribution P(s|0; n). 

4.1. Second and fourth moments of the distribution P{s\0;n) 

Although we do not use them directly in the evaluation of Pn, it is instructive to 
explicitly consider second and the fourth moments of the distribution function P{s\0] n) 
because they provide some insight into some aspects of optimal behavior. 

Consider the second moment of the distribution. One expects that at long times 
the mean square displacement of a particle will be diffusive, that is, that the second 
moment s^(n) for a particle starting from the origin should grow proportionally with 
time n. The overline denotes an average over many trajectories. We differentiate $(/c, z) 
of Eq. ()4.3p twice with respect to k and set = to find that the generating function of 
the second moment obeys 

-Tr\ ^ n-Tr\ az + {1 - a)z^L'^ , 

s^[z) = yz^-s^ln) = ^ 2- (4-6) 

^ {l~az-{l-a)z'^) 

As ^ — >■ 1~ the leading contribution to the generating function is of order 0[{1 — z)^'^], 
which leads to linear growth in time, s'^{n) = 2Dn, with diffusion coefficient 

D = (4.7) 

2{a + {l-a)Tf ^ ' 

Note that for a = 1 and a = this reduces respectively to the standard results D = 1/2 
and D = 1/2/2. 
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One must now distinguish between two situations, the cases T = 1 and T > 1. In 
the first case, when the flight over distance L requires one unit of time precisely as does 
a nearest-neighbour step, Z) is a monotonically decreasing function of a. It is largest 
when a = 0, when all jumps are flights over a distance L requiring a unit of time, and 
achieves its minimal value 1/2 when a = 1, when there are no long-range flights. This is, 
of course, the behaviour one would expect since in the case T = 1 the particle is always 
on the lattice. On the other hand, when T > 1 and the particle spends some portion of 
its time off-lattice, the situation is completely different. Here, D in equation ()4.7|1 shows 
a non-monotonic behaviour as a function of a. Differentiating Eq. ()4.7|) with respect to 
a, we flnd that D has a maximum when 

(r-i)(L2-i)' 

where the superscript "2" stresses that this value of a only corresponds to the maximum 
of the second moment and not necessarily of other moments. Note that < a*-^^ < 1, 
which implies that here, contrary to the case T = 1, the greatest mean square 
displacement is achieved by an intermittent random walk, that is, one that includes 
nearest-neighbour steps as well as off-lattice flights over a distance L. The value of D 
corresponding to a = a^"^^ is given by 

For large L, Dmax may be made very large since Dmax ~ L'^ /T = LV. On the other 
hand, there is a penalty to pay, since s'^{n) = 2Dn is an asymptotic result which holds 
only for sufficiently long times that increases with increasing L. By increasing L and 
optimizing a we can thus make D arbitrarily large, but to observe this regime we will 
have to wait a progressively longer time. 

Consider next the fourth moment s'*(?t,) of the distribution function. Its generating 
function obeys 

TP ^ f ..j^ ^ « + - ">^"'^y , °- + (' - °)^'-" ^ . (4,0) 

ir'o {l-az-{l-a)zTf {1 - az - {1 - a)zTf ^ ' 

The leading behavior of the generating function near z — > 1~ is of order 0[{1 — z)~^], 
so that at sufficiently long times s^(n) is quadratic in n, 

3(a+ 



^ ' {a + {l-a)Tf ^ ' 

We point out again that the coefficient of shows a completely different behavior 
depending on whether T = 1 or T > 1. In the former case, when particles spend all of 
their time on the lattice, it is a monotonic function of a with a maximum at a = 0. On 
the other hand, for T > 1 the coefficient is again a non-monotonic function of a. The 
maximum is achieved at 

^ 2T + LHT-i) 
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For T > 2, the extremum a^^^ is also well within the intermittent regime. 

Note that a*^^-* 7^ a^'^\ (or more precisely, a^^^ < a*^^-*). This alerts us to the fact that 
there is not a single "optimal choice" of a. The choice clearly depends on the property 
or strategy that one wishes to optimize, and for more complex properties it is expected 
to depend on the maximal search time N. We are clear on the particular strategy that 
we wish to optimize, namely, the probability that the target will be detected by at 
least one of the searchers by time N. It is unlikely that either a*^^-* or a^^^ provide this 
optimization. Indeed, we find below that the optimal detection strategy is achieved with 
an A^-dependent value of a that also depends on the jump velocity. 



4-2. Probability of return to the origin and expected number of distinct sites visited 

We now return to the calculation at the end of section 0] to find the generating function 
of the expected value of the number of distinct sites visited by an individual searcher. 
With the substitution of equation ()4.5|1 into equation ()H.12j) . this generating function is 
given by 

^^^^ " {1 - z) {1 - az ~ {1 - a)zT) P{0\0; zY ^^'^^^ 

We thus turn to the analysis of -P(0|0; z). 

We start by noting that equation ()4.4p can be written as 

POO 

P{s\0;z)= / rfte-*J,(Ci,a), (4.14) 
Jo 

where Is{Ci, Cl) is a particular two- variable case of a multi- variable generalized modified 
Bessel function [19], 



/^(Ci,C2, ■••) = - / dkcos{ks)exp 



^ Cm cos{mk) 



m=l 



(4.15) 



in which d = azt, (l = {1 — a)z'^t, and all other (rn are zero. Furthermore, /o(Ci!Cl) 
obeys 

00 

/o(Ci,a) = Io{azt)Ioi{l - a)z^t) + 2 ^ J,i(a^t)/K(l - a)z^t). (4.16) 

1=1 

Subsitituting this latter expansion into the integral on the right-hand-side of 
equation ()4.14|1 with s = 0, performing the integration over t, and using the integral 
representation 

4(p) = ^—-L I e-v^ T„(a;), (4.17) 



X 



where the Tn[x) are Tchebychev polynomials of the first kind, we find that occupation 
generating function for the site s = reads 

P(o|o .) = i /" r + 
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Next, expanding the kernel in a series of Tchebychev polynomials, 

1 1 



1 + azx + (1 - a)z'^y azy/r'^ - 1 



x] 



n=l 



where 



az 



we obtain a closed-form expression for P(0|0,2;) convenient for further analysis. 



P(0|0;z) 



dy 



1-e 



2L 



(4.19) 



(4.20) 



(4.21) 



vraz y_i yr^^T^^ 1 + 2ey + ' 

Note that it is straightforward to perform the integral in equation 1)4.211) in the two 
"pure" limits a = 1 and a = 0. In the first case we have r = l/2;, ^ = (1 — Vl — z^)/z, 
and 



^ Toiy)dy 



l + 25^(-e^)'T,(|/) 



1=1 



tta/I — z^ J -I a/1 — 
On the other hand, when a = 0, we have r = oo, = 0, and 

-/' 



a=0\ 



10; z) 



dy 1 ^ 

1-^/2 1 + Vl-^ 



2T 



(4.22) 



(4.23) 



We notice that for T = 1 the two latter expressions coincide. One might therefore expect 
that P(0|0; z) and the expected number of distinct sites visited will be non-monotonic 
functions of a. 

In the next two sections we explicitly calculate P(0|0; z) and ultimately the desired 
target detection probability for intermittent walks. 



5. Particular case: Lindenberg-Shuler intermittent random walks with L=2 
and T =1,2 

To highlight the optimization procedure, consider first the simple case of intermittent 
random walks with nearest-neighbour steps and jumps over a distance L = 2, first 
studied by Lindenberg and Shuler in [20] and [21]. In this case the integrals in 
equations ()4.4p and (|4.2H) can be performed explicitly. For L = 2 (and T = 1 or 
T = 2), the generating function for the probability of being at the origin is given by 



P(0|0;^; 
where 
Ai= - 

Ao = 



A, 



Ai 



[1 + (1 - a)zT] \ A2 - Ai v^r^ A2 - Ai ^/r^ 



(5.1) 



az 



2[l + {l-a)z^] 



az 



2[1 + (1 -a)z^ 



^ ^ / ^ 8(1 - a)z^ [1 + (1 - a)z^] 



1- 1 + 



\-aW\\^(\-a)z 



a^z"^ 




(5.2) 
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Consequently, in this particular case the generating function of the expected number of 
distinct sites visited is given by the closed-form expression 



S{z) 



(1 + (1 - a)z^) 



(1 - z){l-az-{l- a)zT) 



- Ai) yr^Af yr^A 



.(5.3) 



To find the large-ra behavior of the distinct number of sites visited we need to 
consider the limit z ^ 1^ of the generating function. Recall that we expect a leading 
behavior of the form Sn ~ f{a)n^^'^ because ours is basically a one-dimensional random 
walk, albeit intermittent. Not only do we still need to determine /(a), but we also 
need to ascertain the importance of subsequent terms. A straightforward but tedious 
calculation leads to the expansion 

S(,)-( 2(4-3a) 1 Ajl-a) 1 / 1 \ 

\a + il-a)Tj (1-^)3/2 ay^{4:-3a)y^l-z) y^i-zY^y'^ ' 

Clearly, for any given value of z no matter how close to unity, this "blind" expansion 
in powers of (1 — z) only makes sense if a > since the second term diverges at a = 0. 
This is thus only a valid expansion from which asymptotic behavior can be extracted if 
0<a<l. Ata = one must go back to the original function S{z), set a = 0, and 
then expand: 

2\"' 1 



That a = is a special value that can lead to discontinuities was already noted in [20]. 
It is the only value of the intermittency parameter that absolutely excludes every other 
site from the random walk. 

Inverting S{z) in equations ()5.4|1 and ()5.5j) at times n = N we find that for the 
intermittent random walk with nearest-neighbour steps and jumps to next-nearest- 
neighbour sites the leading large- iV behavior of Sn is 

Sn = f{a)N'" - c{a) + O { ^\ . (5.6) 



Refiecting the discontinuity in a described above, the functions /(a) and c(a) are 
discontinuous at a = 0, 

11/ / I for < a < 1, , , 

f{a)={ \7T{a + {l-a)T)J " (5.7) 




for a = 0, 



and 



4(1 - a) 



for < a < 1, 



cia)={ «i/2(4-3«)V2 - ' (58) 

for a = 0. 

Note that this "discontinuity" should be viewed with appropriate caution, since Sn at 
any fixed N should be a smooth function of a. The discontinuous forms of equations ()5.7j) 
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and (|5.8|) arise because we are describing an asymptotic behavior that is strictly vahd 
only in the limit N —>■ oo. In fact, the use of the a > results is only appropriate as long 
as the second term in Sn does not dominate the first. In any case, these considerations 
require that we recognize the a = discontinuity explicitly. 

The question of interest is whether Sn has a maximum and consequently Pn a 
minimum as a function of a, that is, whether it is possible to design an optimal search 
strategy through the choice of the intermittency parameter. In figure I2ta) we plot the 
asymptotic result for the function Sn, equation ()5.fi|l . and the results of Monte Carlo 
simulations, for = 10^ and T = 1 and 2 as a function of a. We see that indeed 
is a non-monotonic function of a, and that the asymptotic form in equation ()5.(jj) is in 
excellent agreement with the numerical data. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 



a a 

Figure 2. (a) Plot of Sn/VN versus a at N ^ 1000 for T = 1 and T = 2. Solid lines 
are analytic results, equation (|5.6|l . while symbols denote the results of Monte Carlo 
simulations, (b) Plot of the inverse of the non-detection probability, 1/Pn, versus a 
for p = 0.1 and N = 1000. 



Our aim is to analytically determine the value of a which maximizes S^. Since 
we expect a > 0, we use the appropriate form of the coefficients /(a) and c(a) in 
equation ()5.6p . While the coefficient /(a) of the leading term in is non-monotonic, it 
jumps from its smallest value at a = to its largest value at a = 0+. For such small a 
me must retain the second term. Differentiating Sn (both contributions) with respect 
to a to find the maximum readily leads to the optimal value 

^oAN) ^ T (^j^^^y (5.9) 

with additional contributions of lower order in A^. Consequently, the maximal expected 
number of sites visited by the intermittent random walk is attained at a = aopt{N) 
defined by equation ()5.9p . 

max„S„ = (1^) - (4 - T)'/-' (r-/» + -2_) (^) + „ (N"') . (5.10) 
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Let us recap the ingredients of this result and analyze its significance. Because 
the optimal a for the most efficient search is positive definite for any given N, we 
needed to use the asymptotic results for Sn that included the c(a) contribution. The 
actual asymptotic number of distinct sites visited in this optimal search, maXaSiy = 
{32N/7rT)^^'^ , is simply the value of the first term in equation ()5.(j|l with the a > form 
in equation ()5.7|) . when a is set to zero. However, the only reason we could use this 
form is precisely because a is not identically zero. The actual contribution of the second 
term to Sn turns out to be negligible, since it contributes in equation only to order 
N^^^, but we need the existence of this term to obtain the correct asymptotic result. 

In this particular case of a search for an immobile hidden target by an intermittent 
random walk with nearest-neighbour and next-nearest-neighbour steps the best strategy 
can be summarized as follows. If we are ready to wait a sufficiently long time but 
want to be sure that at that moment of time we will attain the lowest possible non- 
detection probability P/v, we choose a = aopt{N) = r(27r/(4 — T)^A^)^/^. Note also 
that this strategy is quite efficient considering that nearest-neighbour and next-nearest- 
neighbour steps are not very different from one other. The enhancement of the expected 
number of distinct sites visited through this strategy is a factor of 2/a/T. This effect 
is exponentiated in the non-detection probability, so that it is dramatically apparent in 
figure m^b). 

6. Optimal search strategy for arbitrary L 

The analysis presented in the previous section for a particular case can be generalized to 
arbitrary values of L and T by following essentially the same reasoning. The algebraic 
derivations that lead to the results are longer and more tedious, but the results are 
similar in form to those obtained previously. The outcome, as we shall see, is that the 
detection probability enhancement can be made much greater than that obtained under 
the restrictions that L = 2 and T = 1 or 2. In fact, it can be improved by orders of 
magnitude. 

As before, to calculate the generating function for the distinct number of sites 
visited we need to determine the generating function for the probability of being at the 
origin. The integral in equation ()4.21|) can no longer be done in closed form for arbitrary 
L and T, but it is possible to carry out an expansion in powers of {1 — z) and retain the 
terms necessary for the computation of Sj\f. As before, the cases < a < 1 and a = 
need to be treated separately because the latter is still a special case. For a = the 
expansion is fairly straightforward and leads to 



The calculation is more elaborate for < a < 1. It is essential to retain not only 
the leading divergent contribution in {1 — z) as z 1~, but also the contribution 
independent of {1 — z). This is because ultimately once again we need the contributions 



Pa=0{0\0-Z) 




1 



1 



(6.1) 
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of 0{N^/'^) as well as 0{N^) to S^- to arrive at the correct optimization. We find 



1 1 



P(0|0;^) = -== -_ + /(a, L) + (v^T^) , (6.2) 

■\/2 (a + (1 - ajT) (a + (1 - ajL^j \/\- z 



where 



TTO; 



and Ti and are the values of the functions r and ^ of equation ()4.2U|) at z 
Substitution into equation ()4.13|) leads to 



S{z) 



2(a + (1 1 21{a, L) (a + {1 - a)L'^ 




(1 - zf (1 - z) 

+ 0\ ,^ \,,A . (6.4) 



which implies that at n = ^ 1 the leading asymptotic behavior of the expected 
number of distinct sites visited by an intermittent random walk is 

:a + (l-«)L2)\'/' /SAT^'/' 



Finally, the integral /(a, L) is not available in closed form, but its important contribution 
to the problem can be estimated. For sufficiently large it is small everywhere except 
for the neighbourhood of a = 0. To extract the leading contribution for small a we 
change the variable of integration in equation (|6.3|) from y to the auxilliary variable 
u = (1 + (1 — a)y)/a, in terms of which 
1-a /-(^-aVa 

/(a, L) 



T^Va Ji {u-l)^{u + l){2-a{u + l)) 

V(l-a)(l + e?^) + 2(aM-l)ef' V^7^(a + (1 - a)L2 ' • 



The singularity at the lower limit of integration m = 1 is integrable for any value of a. 
Furthermore, the integrand vanishes sufficiently fast to insure that it is integrable as 
u ^ oo. Hence, we may set a = everywhere in the integrand and we can extend the 
upper integration limit to infinity. This provides the leading small-a contribution, 

L) ~ g{L) = r f coth(L0) - J coth(0)') . (6.7) 

TTy/a Jo smh(0) \ L J 

The subsequent reasoning now proceeds exactly as in the previous section. We wish 
to determine the value of a that maximizes S^- Again, we expect that this a is small 
but strictly positive, so we focus on the result ()6.5|) . Maximizing with respect to a 
we find that 

- Vl/3(^y _ 1)2/3 (87rAr)V3' ^^-^^ 
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The resulting asymptotic form for the distinct number of sites visited finally is 



maXa^AT = (LV) 



1/2 



TT 



1/2 



iN 



71 



1/6 



where 



(471 + 1) 
27l2/3 



Vg\L)L\LV 



(6.9) 



(6.10) 



This is the generalization of equation ()5.10|) . It is valid when aopt{N) <^ 1, that is, 
when ^ g'^{L)L^, which insures that the first term in max„5'Ar is dominant. Note 
that the leading contribution to Sn is again the outcome of setting a = in the first 
term of (|6.5|) . The second term ultimately does not contribute to leading order in A^, 
but, to stress once again, we were required the use the a > result ()6.5|) because aopt{N) 
is positive definite. 

In figure a) we plot our analytic result ()6.5|1 using the estimate ()6.7j) (lines) as well 
as the results of Monte Carlo simulations (symbols) as a function of a. The results are 
shown for L = 5 and various values of and T, and the agreement is again excellent. 
The larger symbols in the figure denote the positions of the maxima, which are well 
captured by our analytical curves. Figure Efb) shows the inverse of the associated 
non-detection probability for various values of N. 
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Figure 3. (a) Plot of Sn/\/N versus a for L = 5 at iV = 10000, for T = 1, 2, 3, 4, 5. 
Solid lines are analytic results, equation (|6.5|l with (|6.7|l . while symbols denote the 
results of Monte Carlo simulations. The larger symbols denote the positions of the 
maxima, (b) Plot of the inverse of the non-detection probability, 1//V, versus a for 
p = 0.01 and N = 8000, 10000, and 12000. 



The essential result of this section is the very large enhancement of the target 
detection probability that is possible by the appropriate iV-dependent choice of the 
intermittency parameter compared to the outcome for the "pure" random walks with 
a = or a = 1. The leading term in equation ()6.9|) contains the factor V LV which 
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may be large and may therefore lead to a substantial decrease in the probability that 
the target has not been detected up to time N. For example, for L = 5, T = 1, and 
= 10000, with a density of searchers as low as p = 0.01 (cf. figure Efb)), the non- 
detection probability for the "pure" cases is P/v ~ 0.2. On the other hand, if we choose, 
with a = aopt(lOOOO) ~ 0.07, the probability that the target has not been yet detected 
up to this time is P/v ~ 0.0003, i.e., three orders of magnitude smaller. Note finally 
that in one dimension for fixed L the optimal search strategy involves a progressively 
smaller fraction of nearest neighbour steps as is increased. 

7. Conclusions 

We have studied the search kinetics of a single target hidden on a one-dimensional 
regular lattice by a concentration of "searchers" performing intermittent random walks; 
that is, at each tick of the clock, each searcher has a choice of stepping on a nearest 
neighbour with probability a, or flying off-lattice in either direction with velocity V 
over a distance L with probability (1 — a). We have determined the probability P/v that 
the target remains undetected up to the maximal search time A^, have established that 
P/v is a non-monotonic function of a, and that in fact it has a sharp minimum at an 
A^-dependent value of a. Consequently, we have shown that the search efficiency can be 
dramatically enhanced, by orders of magnitude, by choosing a appropriately. 

We close this paper with some observations and problems for future research. If we 
are at liberty to choose L, we may increase the number of distinct sites visited within 
a given time interval A^ by picking an A^-dependent value of L, cf. equation ()(i.9|) . On 
the other hand, equation ()6.9|) is a parabolic function of L, which implies that there 
should be an optimal L for any fixed A^. To wit, making L too large for fixed A^ would 
be counter-productive since the searchers would spend most of their time fiying off- 
lattice. On the other hand, a value of L that is too small would detain the searchers in 
already well-explored regions and delay the exploration of new ones. This issue requires 
a more subtle analysis and a more intricate optimization procedure, to be discussed 
elsewhere [15]. 

Our model can be straightforwardly extended to higher dimensions. In two 
dimensions there is also a value of a that maximizes and that depends on the 
maximal search time A^, as seen in the numerical results shown in figure IH In three 
dimensions and higher, we expect an A^-independent value of a to maximize 5'/v. 

The results in this paper were obtained for the case of perfect detection (g = 1), that 
is, once a searcher lands at the site occupied by the target it detects it with certainty. 
It is easy to show that in the general case when detection of the target upon encouter 
is less than certain (0 < g < 1), the probability that the target remains undetected up 
to time A^ is 




(7.1) 
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Figure 4. Results of Monte Carlo simulations on a two-dimensional square lattice. 
Plot of SNln{N)/N versus a for different A''. Symbols denote the results of the 
simulations, and the solid lines are just guides for the eye. 



where Qiv is determined via its generating function, 

oo 

«w = = — |,p(oio;:).i-,i i:^Mo;^)- (^-2) 

In low dimensional systems, i.e., in one and two dimensions, -P(0|0; z) diverges as 2; ^ 1~, 
and consequently, in the limit — > 00, Qj^ converges to so that the leading large- 
behavior will be independent of q provided that q > 0. In three and higher dimensions 
the leading behavior will depend on q. 

An interesting extension of our model may be a situation in which the target itself 
moves randomly. It is known, however, that if this motion is diffusive (or sub-diffusive) 
in low dimensions, the long-time asymptotic form of the probability Pjv is exactly the 
same as when the target is immobile (see [22] and also [23,24]). Consequently, in low 
dimensions in the largc-iV limit, the search for a diffusive (or sub-diffusive) target by 
searchers performing intermittent random walks will proceed in exactly the same way 
as determined in this paper for the case of an immobile target. 

Finally, we remark that a sensible extension of this model would be to introduce 
the possibility of jumping, with a given velocity over arbitrary distances / = 1, 2, 3, . . ., 
with probability pi. A robust approach would be to search for the distribution pi that 
minimizes the non-detection probability P^- 
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